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ABSTRACT 

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum 
of the Ultra-compact HII region G29. 96-0.02. Spectral lines from species including 13 CO, CO, [CI], and [Nil] are detected. A sparse map of the 
[Nil] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS 
and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. 
The quality of the spectrum longwards of 100 yum allows us to fit a single temperature greybody with temperature 80.3 + 0.6 K and dust emissivity 
index 1.73 + 0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 M G for the clump containing the HII region. 
The clump's bolometeric luminosity of 4 x 10 6 L G is comparable to, or slightly greater than, the known O-star powering the UCHII region. 
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1. Introduction 

Ultra-compact HII (UCHII) regions are small nebulae that sur- 
round massive, young stars that are still emb edded within a natal 
cloud dChurchwel]|2002tlHoare et alJ2007t) . UCHII regions rep- 
resent an important evolutionary stage in the formation of mas- 
sive stars. They may correspond to the still embedded accretion 
phase while the star is already ionising its surroundings. Their 
study is critical to our understanding of how stars manage to 
reach high masses despite the high pressure generated by their 
growing HII regions. UCHII regions are amongst the brightest 
Galactic objects at submillimetre wavelengths and can be used 
as a probe of star formation in distant galaxies. Thus, they make 
excellent candidates with which to probe all phases of massive 

star formation and with which to test new instruments; 

G29.96-0.02 is an UCHII region (IWood & Churchwelll 



ll989HPeetersetalJ|2002l) located at a distance of 8.9+°^ kpc 



(ISewilo et al.1 120041) . The bright I R source at the centre 
of G29.96-0.02 is an Q5-6 star (Wats on & Hanson! Il997t 
iMartfn-Hernandez et alj|2003l). with a lumin osity of Lboi ~ 3 — 
4 x 1 6 L Q for the[ Watson & Hanso3 (1 19971) temperature limits 
at the ISewilo et al.l (120041) distance. However, models indicate 
that a cluster of young stars must be present in G29.96-0.02 to 
account for its observed luminosity (Lumsden et al. 2003). 

At the head of the cometary shaped radio emission 
(IWood & C hurchwelll 119891) is a small cluster of hot molec- 
ular cores (~300K) with masses of 3-1 1M (Beuth er et alj 
2007). These arcsec-sized structures are embedded within an 
arcmin-sized submillimetre clump called G29. 956-0. 017SMM 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



(iThompson et al.l 120061) . iFev et all (1 19951) suggested that the 
cometary shape of the UCHII region is caused by its expan- 
sion westwards into the denser part of the surrounding clump. 
This clump is part of a larger complex that has been mapped by 
ifersc/ie/- SPIRE/PACS as part of the Hi-GAL key programme 
dMolinari et alJ r2010, Beltran et al. in prep.). 

In this paper we present new submillimetre spectra taken 
towards G29.956-0.017SMM (hereafter G29) with the ESA 
Herschel Space Observatory (Pilbrattet al. 2010) as part of the 
"Evolution of Interstellar Dust" key programme (Aberge l et al.l 
2010). The data were taken with the SPIRE Fourier- Transform 
Spectrometer (FTS). The SPIRE instrument, its in-orbit per- 
formance, and its scientific capabilities are described by 
Grif fin et aD d2010h . and the SPIRE astr onomical calibrat ion 
methods and accuracy are outlined bv lSwinvard et alj d2010l) . 

2. Observations 

G29 was observed with the high-resolution mode of the SPIRE 
FTS on 13 September, 2009 at 19:49 (Herschel observation 
ID, 1342183824). Two scan repetitions were observed giving 
an on-source integration time of 266.4 seconds. The point- 
ing centre was at a Right Ascension and Declination (J2000) 
of 18 A 46 m 04.07\ -02°39'21'.'88 and Galactic coordinates of 
29.9561, -0.01726. The unapodized spectral resolution was 
1.2 GHz (0.04cm~'), after apodization (using extended Norton- 
Beer function 1 .5, iNavlor &T ahic 2007) this became 2. 17 GHz. 
The SPIRE FTS measures the Fourier transform of the 
source spectrum across short (SSW, 194-313yum) and long 
(SLW, 303-671 /mi) wavelength bands simultaneously. Each 
waveband is imaged with a hexagonal bolometer array with pixel 
spacing of approximately twice the beam-width. The FWHM 
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Table 1. Best Fit Line Parameters 
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Fig.l. SPIRE FTS SLWC3 (red) and SSWD4 (blue/grey) 
spectra of G29. Upper: On-source spectrum with continuum, 
point-source calibration. Middle: As above with continuum sub- 
tracted. Major lines are annotated. Lower: Mean off-source 
spectrum, extended calibration. See text for details. 



beam-widths of the SSW and SLW arrays vary between 17- 
21" and 29-42" respectively. The source spectrum, including the 
continuum, is restored by taking the inverse transform of the ob- 
served inte rferogram. For deta ils o f the SPIRE FTS and it s cali- 
bration see lGriffmet~aT1 (l2010) and lSwinvard et all (120101) . 



3. Results 

The upper panel of Fig. Q] shows the apodized SLW (red) and 
SSW (grey) spectra. Only the data from the central bolometer 
(C3 for the SLW and D4 for the SSW) as calibrated for a point 
source are shown for each spectrum. The spectra are dominated 
by the thermal continuum. Superimposed on this are a series 
of bright lines, the most noticeable of which are the ladder of 
CO lines. The shape of the continuum was estimated by mask- 
ing the CO lines and performing a linear regression of the form 
log S v = C + p log v to both spectra. The data were best fit by 
power-laws with indices of pssw = 2.71 and pslw - 3.40 re- 
spectively. There is a disconnect between the two spectra. The 
offset in power-law constants was Cssw - Cslw - 0.164, equiv- 
alent to a linear factor of 1.45. The SSW spectrum was shifted 
upwards by this margin and replotted as the blue spectrum under 
the assumption that the discontinuity was due to differing struc- 
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Fig. 2. Population diagram for 12 CO and 13 CO towards G29. 



ture within the SLW and SSW beams. The power-laws for each 
spectrum are plotted as black lines. 

The Leve nberg-Marquardt least-squares minimisation pack- 
age MPFIT (iMarkwardtl 120091) was used to simultaneously fit 
a catalogue of lines and an 8th-order polynomial to each of 
the SLW and SSW spectra. It was assumed that the lines were 
Gaussian and that the linewidths were not resolved. The mid- 
dle panel of Fig. [TJshows the spectra after the polynomial back- 
ground has been subtracted. The best-fit line shapes are shown 
in red and blue for the SLW and SSW bands respectively. These 
show the same ladder of CO lines as the upper panel as well 
as their 13 CO count erparts. The positio n of the I2 CO and 13 CO 
lines are annotated (Pickett et al. 1998). The brightest non-CO 
lines are the 492 and 809 GHz lines of [C i]. The 809 GHz [C i] 
line is blended with the 806 GHz J=7-6 line of 12 CO Table [T] 
lists the transitions, reference frequencies, frequency offsets, and 
peak flux densities for each significant line fit. 

Figure|2]shows the rotational po pulation diagram for | 2 CO in 
green and 13 CO in blue. Following Goldsmith & Langer ( 1999) 
we fit a rotational temperature to each species. The weighted lin- 
ear fits and results are plotted for each species. The quoted errors 
are the errors on the fit and do not include uncertainties in beam 
size, calibration, and sub-structure within the beam. The 13 CO 
lines come from a greater depth of material than the 12 CO lines 
and will be more indicative of the temperature of the interior 
of the clump. The 7-6 transition of 13 CO and the 9-8 transition 
of I2 CO appear to be lower than the trends. This could be due 
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Fig. 3. Left: False-colour map of G29 showing GLIMPSE 4.5 
(blue) and 8-//m (green) (Benjamin et al. 2003) and MAGPIS 
20cm (red) dHelfand et aUl200 6). Right: [N n] line flux towards 
G29. The circles show the position of the bolometers on the sky. 
The space between them has been interpolated using inverse 
distance weighting. The X m arks the position of the G29.96- 
0.02 radio source dWood & Churchwelll [19891 18''46'"03.96 s , 
-02°39'21'.'5 ). The diamond marks the position of the Nil peak. 
The contours are SCUBA 850/mi emission (5, 10, 20, 50 and 
100 times the local off-source rms) showing the extent of the 
G29.956-0.017SMM clump. 



to the structure in the beam changing between transitions (the 
beam FWHM changes by a factor of 2.5 between the CO 4-3 
and 13-12 transitions). 

The region around G29 is shown in the left panel of Fig. 
[3] as a false -colour image - blue and green are GLIMPSE 4.5 
and 8-um dBeniamin etaf] 12003) and red is MAGPIS 20cm 
dHelfand etal.ll2006l) - and archival SCUBA 850/mi contours. 
The X marks the positio n of the VLA radio e mission associated 
with the UCHII region dWood & Churchwelllll989l) . The spec- 
tra measured towards positions surrounding G29 show several 
notable differences to the on-source spectra. The lower panel of 
Fig. [TJ shows the mean background-subtracted spectrum for co- 
aligned off-source bolometers. The 835 GHz line of CH + is seen 
in absorpt ion in the SLW b and and appears to be probing the 
ISM (see Navl or et al.l2010l for a detailed decomposition of the 
CH + line for this source). The deepest absorption feature is coin- 
cident with the 1.232THz line of HF (Neufeld et al. 2010). The 
higher-order I2 CO and 13 CO lines are not detected in the SSW 
band, but the 1.46 THz fine structure line of [N n] is strongly 
detected. This line is present in the extended-source calibration, 
but not in the on-source point-source calibration. 

A fit to the [N n] line was performed for all bolometers of 
the SLW array. The results are plotted on the right panel of Fig. 
[3] The coloured circle markers show the bolometer positions and 
the measured line intensities at those positions. For these early 
results the bolometer array was not offset to create a fully sam- 
pled map. To compensate for this the pixels between the bolome- 
ters have been interpolated using an inverse distance weighting 
(modified Shepard's) method. The southern [N n] peak is coin- 
cident with a region of 20cm emission that is bounded to the 
north by a 8 -/mi filament suggesting that it is a separate HII 
region to the G29 UCHII region. Its diameter is ~V (2.5 pc at 
8.9 kpc). The morphology around G29 is clearly complex and a 
more detailed, fully-sampled study will be required to disentan- 
gle the various components. As stated, this preliminary map was 
not fully sampled, but it does shows the potential of the FTS as 
a mapping spectrometer. 



4. Spectral energy distribution 

Figure [4] shows the spectral energy distribution (SED) measured 
towards G29. SPIRE FTS SLW-C3 and SSW-D4 spectra are 
shown in the same colours (red and blue) as Fig. [TJ A search 
was made for archival data coincident with G29. For consistency 
with the FTS data the downloaded spectra and data points were 
converted into intensity units by dividing by their wavelength- 
dependent beam area or, where relevant, their aperture area. The 
archival data is plotted in Fig.|4]and described in the Figure cap- 
tion. 

Figure [4] does not include the relative correction applied in 
Fig. [TJ It was found that the FTS SLW data were marginally 
higher that the SCUBA 450^m data point and that the first point 
of the FTS SSW spectrum was significantly higher than the ISO 
Longwave Spectrograph (LWS) final data point. It was assumed 
that the three sections of the long wavelength SED (ISO LWS, 
FTS SSW, FTS SLW) followed the same single temperature 
greybody and that the offsets between them were due to differ- 
ences in their absolute calibration. A greybody equation was fit- 
ted to these three spectra and offsets between them accounted for 
by giving the flux from each section a multiplier coefficient. The 
greybody had the form 



I v = C x B v (T){l-e- {vKf ) 



(1) 



where C x is the coefficient, B V (T) is the Planck function at tem- 
perature T, v c is the frequency where the emission becomes op- 
tically thin, and/? is the dust emissivity index. The data were re- 
sampled into evenly spaced bins in log-wavelength space to pre- 
vent the differing density of data points between the ISO LWS 
and the SPIRE FTS corrupting the minimisation. 

The best-fit to the SED was T = 80.3 ±0.6 K, B = 1 .73 ±0.02, 
and v c = 20.0 ± 0.5 THz (equivalent to A = 15.0 ± 0.4 /mi). 
We note that the SED temperature is similar to the 13 CO rota- 
tional temperature from Fig. [2] The multiplicative coefficients 
for the FTS SSW and FTS SLW were found to be 9.3 ± 0.08 
and 4.04 ± 0.06 respectively, the ISO LWS coefficient was fixed 
at 1.0. The SPIRE FTS spectra have been aligned to the ISO 
LWS calibration. A single temperature greybody well fits the 
data longwards of ~ 40 //m. However, fitting a single compo- 
nent greybody to a complicated source like G29 can only yield a 
broadly characteristic temperature. The relatively flat SED short- 
wards of the peak shows that there must be several hotter temper- 
ature components. Nevertheless, Fig.|4]does show that the long 
wavelength the slope of the SPIRE FTS data is consistent with 
the slope of the literature data. 

In the following paragraphs we adopt the distance of 8.9 kpc 
(Sew ilo et al.l 120041) for G29 and assume that the majority of 
the emission comes from a region comparable to, or smaller 
than, the Hershel beam at 250 /mi (0=18") - (at 850 /mi, 60% 
of the extended emission towar ds G29 is within the central 
SCUBA 14.7" beam, (Thompson et al.ll2006l) ). Based on the fit- 
ted dust temperat ure, G29's published 850/mi peak flux density 
(Thompso n et al.l 12006). and making typical assumptions (e.g. 
Kirk et al. 120051) we estimate the mass of the G29 clump to be 
M ~ 1500 M . As with the fitted temperature, the actual mass 
will depend on the internal profile/geometry of the clump, but 
we note that this mass is similar to the median mass (940 M Q ) of 
the infra-red dark cl ouds (IRDCs) from wh ich high-mass stars 
are believed to form (Rathbor ne et al.l|20 06). 

The luminosity integrated under the fitted greybody in the 
range 2-2000 /mi is Lo usl - 61 D 2 9 2 L Q where D is the distance 
to the source in kpc and 9 is the diameter of the emitting area. 
The above assumptions give a luminosity of Lo„„ = 1.6xl0 6 L o . 
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Fig. 4. The G29 spectral energy distribution. SPIRE SWS-D4/SLW-C3 spectra ar e shown by the blue/re d lines. The light and 
dark grey spectra are archi val data from the ISO Short Wavelength Spectrograph dde Graauw et aLj 19961 SWS) and ISO Long 
Wavelength Spectrograph (Cle22 et al. 1996, LWS) respectively and were originally published bv iPeeters et alJ ( 20021). Archival 
MIPS SED mode data (Rieke et al. 2004) are shown by the green line. The error bars show IRAS 12. 5-100 /im (Helou & Walker 
1988, IRAS 18434-0242 is coincident with G29), SCUBA 450 and 850/im dThompson et alj|2006l) . and IRTF 1.3irun dChini et al. 
1986) data points. The FTS SED data have been shifted to the same calibration as the ISO LWS data. A best fit SED is shown by 
the dashed line. 



This agrees with the luminos ity of 10 6 L Q e stimated from IRAS 
measurements alone (IWood & Churchw ell 1991). Likewise, the 
bolometric luminosity in the range 2-2000 /^m, interpolating to 
the fitted SED at A > 650//m, is L bo , = 154D 2 6 2 L . The 
greybody luminosity is ~40% of Lboi- At the assumed distance 
Lboi = 4.0 x 10 6 L . The Lboi of the clump containing the UCHII 
should equal the luminosity of the driving sources if the region 
is in equilibrium. Our calculated Lboi is on the upper limit of 
the range of luminosities for the identified single O-star in G29 
(3 - 4 x 10 6 L , Watson & Hanson 1997) which may support the 
idea that the luminosity from more than one star is powering the 
reprocessed SED dLumsden et al. 2003). 

5. Conclusions 

We have presented new SPIRE FTS 190-670yum spectra of the 
submillimetre clump G29.956-0.017SMM which contains the 
G29.96-0.02 UCHII region. The impressive capabilities of the 
SPIRE FTS have allowed us to simultaneously observe both the 
dust continuum and prominent spectral lines towards G29. We 
have conducted basic line-fitting and shown the distribution of 
[N n] emission towards G29. While this preliminary map was not 
fully sampled it does show at least one other HII region neigh- 
bouring G29 and demonstrates the amazing potential of the FTS 
as a mapping spectrometer. We have reconstructed the SED of 
G29 using the FTS spectra and archival measurements and have 
shown that the FTS calibration is broadly consistent with earlier 
observations. The combined data-set allowed us to fit a precise 
greybody with T — 80 K and /3 - 1.73. Based on a distance 
of 8.9 kpc we estimated the mass of G29 to be approximately 
1500 M . The calculated luminosity of the G29 clump is slightly 
greater than the known O-star at the centre of the UCHII region. 
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